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Constraints on the parameters in the one- and two-loop pion-pion scattering am- 
plitudes of standard chiral perturbation theory are obtained from explicitly crossing- 
symmetric sum rules. These constraints are based on a matching of the chiral ampli- 
tudes and the physical amplitudes at the symmetry point of the Mandelstam plane. 
The integrals over absorptive parts appearing in the sum rules are decomposed into 
crossing-symmetric low- and high-energy components and the chiral parameters are 
finally related to high-energy absorptive parts. A first application uses a simple 
model of these absorptive parts. The sensitivity of the results to the choice of the 
energy separating high and low energies is examined with care. Weak dependence 
on this energy is obtained as long as it stays below ~ 560 MeV. Reliable predictions 
are obtained for three two-loop parameters. 
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I. INTRODUCTION 



A method for the determination of the parameters appearing in the pion-pion scattering 
amplitudes produced by chiral perturbation theory!'! has been proposed in!. We demon- 
strate its practicability in the present work. The method is based on sum rules derived from 
exact analyticity properties combined with crossing symmetry. A similar approach has been 
developed and applied in!. Two features distinguish our method: 

(1) Our constraints on the chiral parameters are obtained by matching the true ampli- 
tudes and their chiral approximations at the symmetry point s = t = u = 4M 2 /3 of the 
(s, t, w)-space (s, t and u are the standard Mandelstam variables). More precisely we equate 
coefficients in the Taylor expansions of both amplitudes at the symmetry point. 

(2) The Taylor coefficients are expressed by means of crossing-symmetric sum rules 
as integrals over absorptive parts. These integrals are decomposed into exactly crossing- 
symmetric low- and high-energy components. If A is the energy separating low and high 
energies, we assume that A can be chosen in such a way that the absorptive parts can be 
approximated by chiral absorptive parts below A and that they are obtained from reliable 
experimental data above A. 

The last assumption does not fix A precisely, whereas the form of the conditions con- 
straining the chiral parameters depends explicitly on A. Our method is consistent only if 
these conditions lead to values of the parameters depending weakly on A. This stability 
of the parameters with respect to a variation of A is one of our main concerns. To settle 
this point we work with a simple model of the true absorptive parts defined for all energies 
and at least qualitatively consistent with the data. The separating energy A can be pushed 
down to the elastic threshold 2M n . We have constructed the constraints for the parameters 
of the one- and two-loop pion-pion amplitudes of standard chiral perturbation theory for 
values of A 2 ranging from 4M 2 to 40M 2 . Our assumptions turn out to be valid as long as 
4M 2 < A 2 < 16M 2 (= (560 MeV) 2 ). 

Before we continue we have to define our notation and the chiral parameters we are 
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using. The s-channel isospin I amplitude T J (s, i, m) has the partial wave expansion 

T^s, t, u) = ]T(2/ + l)tj (s)P, (l + , 7 = 0,1,2. (1.1) 

The pion mass has been set equal to 1, consequently s + t + u = 4. Below the inelastic 
threshold 



ti(8) = J— f &Usm6 l {8), (1.2) 

where the phase shift 5/(s) is real. 

The three amplitudes T 1 are obtained from a single function A(s;t, u) : 

T°(s, t, u) = 3A(s; t, u) + A(t; u, s) + A(u; s, t), 

T\s, t, u) = A(t; u, s) - A(u; s, t), (1.3) 
T 2 (s, t, u) = A(t; u, s) + A(u; s, t). 

In chiral perturbation theory A is the sum of a polynomial A po1 and an analytic function 
-A cut . At 6th order in the momenta and the quark masses, each component of A can be 
written as 

A® = X 2 A 2 + A 4 A 4 + X 6 A 6 , (1.4) 
where A = M n /F n . In standard chiral perturbation theory we have! 

OZ7T 

Al°\s; t,u) = -^- [fei + b 2 s + b 3 s 2 + b 4 (t - uf] , (1.5) 



32tt 

Al°\s;t,u) = -^-\b 5 s 3 + b 6 (t-u] 

OZ7T L 



The coefficients fej are the parameters we want to determine. Their relationship to the 
coupling constants in the 4th- and 6th-order effective chiral Lagrangian is given in!. We 
do not use these relations here, remarking only that b\, b 2 , b% and 64 are sums of 4th- and 
6th-order terms whereas 65 and bg are 6th-order parameters. 
The analytic part of A^ has the form 



(1.6) 

+s[T a (t)/a(0+^a(w)/a(u)]} 

where R Q (s), S a (s) and T a (s) are polynomials containing 4th- and 6th-order terms and 
/ a (s) are analytic functions with a right-hand cut [4, oo). The polynomials are obtained 
from formulae (2) and (3) in!: the functions used there have been relabelled f\(s) = J(s), 



fa(s) = K a -i(s), a = 2, . . . , 5. The coefficients of the polynomials in ( |1.6|) depend linearily 
on 6j. 

Our basic sum rules are obtained from subtracted dispersion relations and contain un- 
known subtraction constants. For this reason the number of useful constraints is smaller 
than the number of parameters, two constraints at 4th order and four constraints at 6th 
order. They are used for the determination of b^ , 64 , bf^ , bf* , bf^ and bf^ at fixed values 
of bi and b% ■ Four parameters fulfil the criterion of near A-independence: 6 4 4 \ 63 , b^ 
and fog 6 '* . In spite of the fact that our input modelling the experimental data is rather crude, 
the values we get for the stable parameters are entirely compatible with the results derived 
ini. 

The construction of our constraints is outlined in the next two Sections. The inputs are 
defined in Section IV and the chiral parameters are evaluated in Section V. The results are 
discussed in Section VI. 



II. SUM RULES FOR TAYLOR COEFFICIENTS AT THE SYMMETRY POINT 

The full amplitudes T 7 (s, t,u) as well as the chiral amplitudes T^(s,t,u) are real and 
regular in the Mandelstam triangle 0<s<4, 0<t<4, 0<m<4. This implies that their 
Taylor expansions at the symmetry point s = t = u = 4/3 converge in a neighbourhood 
of that point. We assume that the 2n-th order chiral amplitudes approximate the full 
amplitudes in that neighbourhood up to O (p 2(n+1 )^ corrections. If this assumption is correct 
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we can fix the parameters appearing in the chiral amplitudes by requiring that T 1 and T£ 
have identical conveniently truncated Taylor expansions at the symmetry point. We adopt 
this strategy and explore its implications. 

To write down Taylor expansions for the T 1 , s, t and u have to be replaced by two inde- 
pendent variables. Generally the series one obtains are constrained by crossing symmetry. 
It is advisable to resort to a procedure which ensures that these constraints are automati- 
cally fulfilled. To this end the three amplitudes T 1 are replaced by three totally symmetric 
amplitudes Gi (i = 0, 1, 2) expressed in terms of two new variables x and y which are totally 
symmetric and homogeneous in s, t and t/i: 



1 / 

x = {st + tu + us) 

16 



y = — stu. 
y 64 



(2.1) 



The Gi have been introduced by Roskiesi and are defined as follows: 
1 



G Q (x,y) 



T°(s, t, u) + T°(u, t, s) + T°(u, s, t) 



t — u 



u 



(2.2) 



s - t 
1 



t — u 



s — u 



+ 



t — u 



T\t,u,s) T l {u,t,s) 



u — s 



t - s 



+ 



u 



T^u^^) T\s,u,t) 



s — t u — t 

The T 1 can be obtained from the Gi (see eq. (2) in!); their symmetry implies crossing 
symmetry of the T 1 . 

The Gi are regular around the image (x s ,y s ) of the symmetry point in the (x,j/)-space 
(x s = —1/3, y s = 1/27) and have convergent Taylor expansions in x and y at that point. 
They are not constrained by crossing symmetry. We equate their first coefficients with the 
corresponding Taylor coefficients of the symmetric chiral amplitudes Gf, the latter being 
given by Q with T 1 replaced by T£. 

The Taylor expansions of the amplitudes G^ obtained from the 0(p 2n ) amplitudes 
T( 2n ) are truncated at an order obtained by power counting from the definitions (|2.2j ). When 
identifying truncated Taylor series at a given order one has to remember that a chiral 



Taylor coefficient gets contributions from all higher orders of the chiral expansion. Therefore 
the conditions one obtains for the 2n-th order amplitudes hold only up to higher order 
corrections. 

For the fourth- and sixth-order amplitudes we get the following constraints: 

Gf\x s , Vs) + 0( P 6 ) = Gi(x s , y s ), i = 0, 1, 2, (2.3a) 
8G^ r)G 

-^( Xs7 y s ) + 0(p 6 ) = -Q^(x.,y.), (2.3b) 



G?\x„y a ) + 0(p B ) = Gi(x a ,y 8 ), i = 0, 1, 2, (2.4a) 
' / x s ,y s ) + 0(p 8 ) = ^(x s ,y s ), i = 0,l, (2.4b) 



OX ox 
^-(x„, y.) + 0(p 8 ) = ^(x s , y s ). (2.4c) 

The number of constraints obtained in this way at 2n-th order is equal to the number of free 
parameters appearing in the 2n-th order chiral amplitudes (4 parameters and 4 equations 
(2.3) at 4th order, 6 parameters and 6 equations (2.4) at 6th order). An obvious drawback of 
conditions (2.3) and (2.4) is that, the symmetry point being unphysical, the right-hand sides 
are not directly measurable quantities. Fortunately, there are dispersion relations to get rid 
of this dimcultyi. These relations are consequences of the exact analyticity properties of the 
Gi as functions of x and y. They can be written in the following way: 

1 1-5,2 



1 f°° 

Gi(x, y) = (1 - d i2 )Gi(xo, y ) + - / da 

7T Ji 



Bi(a,r). (2.5) 



K(<r,x,y) K(a,x ,y ) 
There are once-subtracted dispersion relations for G and G\, the subtraction being per- 
formed at (x ,y ), whereas G 2 obeys an unsubtracted dispersion relation. The integration 
variable a is an energy squared. The denominator function K(a, x, y) has a simple expression 
in terms of the Mandelstam variables: 

K(a,x,y) = {a- s)(a-t)(a - u) = a 2 (a — 4) - 16ax - 64y, (2.6) 

if (s,t,u) is a pre-image of (x,y). The second expression results from the first one and the 
definitions (2.1). The Bi are linear combinations of s-channel absorptive parts A 1 : 
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B Q (s,t) = -(s-t)(2s-4 + t) \A°{s,t) + 2A 2 {s,t) 



BAs, t) = -(3s - 4) \2A°(s, t) - 5A 2 (s, t) 
6 L 



+ 



(s-t)(2s-A + t) 1 
2t - 4 + s 2 



(2t - 4 + s) 



B 2 {s,t) 



1 



2A°(s,t) -5A 2 (s,t) 



3 3s -4 4l . . 

+ A 1 (s,t). 

2 2t - 4 + s V ' ; 



(2.7) 



They have to be evaluated in (2.5) at a a-dependent squared momentum transfer r: 



t 0" 



(a-4) 2 



16 



-(aa((j - 4) - 16(ax - y )) 



(2.8) 



a + 4a 

where a is the slope of the straight line connecting the point (x, y) to the subtraction point 
(x ,y ): 

y-yo 



X — Xq 



(2.9) 



Details on the derivation of (2.5) are given ini. At fixed (xo,?/o)> (2-5) holds true if the 
slope a belongs to a complex neighbourhood of the origin. For our purpose (xq, yo) has 
to be chosen in such a way that (2.5) provides a representation of Gi(x,y) in a suitable 
neighbourhood of the symmetry point (x s , y s ). According toi this is the case if yo = y s and 
-72 < x < 3x s /2. 

If er is large enough, r(er) is real and the point (s = a, t = r(cr)) belongs to the physical 
s-channel. If a is close to 4, r(cr) can be complex but the point (s = cr,t = t((j)) is inside 
the large Lehmann ellipse. This implies that the absorptive part A 1 (a, r) is either a physical 
quantity or is obtained from physical partial wave absorptive parts through a convergent 
partial waves expansion. Consequently the dispersion integral in (2.5) is itself a physical 
quantity and this relation produces a representation of G 2 (x,y) in terms of measurable ab- 
sorptive parts. This does not hold for Go(x,y) and Gi(x,y) because their values at the 
subtraction point which enter into (2.5) are unphysical with our choice of (xq, yo). However, 
by computing derivatives at fixed (x ,y ), (2.5) gives expressions for dGi/dx and dGi/dy 
at (x s ,y s ) which do not involve Gi(xo,yo)- Therefore the dispersion relations provide ex- 
pressions in terms of observable absorptive parts for the right-hand sides of eqs (2.3a, i=2), 
(2.3b), (2.4a, i=2), (2.4b) and (2.4c). 



If one wants to keep eqs (2.3a) and (2.3b) for i = 0,1, one can use a dispersion relation 
connecting Gi(xo,yo) to the value of G{ at the elastic threshold x = y = (image of 
s — 4, t — u — 0). This value is determined by the S-wave scattering lengths, and one 
obtains expressions for Gi(x s ,y s ), i = 0,1, in terms of physical quantities. However, as the 
experimental scattering lengths are poorly known at present, this would not lead to stringent 
constraints on the chiral parameters. Therefore we restrict our discussion to the constraints 
involving only dispersion integrals over absorptive parts. 

We end this Section with the expressions of the derivatives dGi/dx, i = 0,1, and dGo/dy 
resulting from (2.5): 

dd 32 r°° q 2 (cr-2) + 32y 8 1 

{x s ,ys) = — da—— — — - — — — -Bi(a,T), i = 0,l, 



&r VA '^ 77 7 4 * (K(a,x s ,y s )y {a - r) (2a - 4 + r) 
8G 0/ \ 64 roc 1 f pffo _ g) _ iq x i 

2(o- 2 (q-2) + 64 j/s ) 1 \ 1 

" a^[a{a - 4) 2 - 256yJ 1 /a 3 \W A + 2A >)^ T )j- 

In these formulae, r is given by (2.8) at a = 0. As a consequence, the integrands do not 
depend on xq, the position of the subtraction point on the line y = y s . This is an advantage 
of our choice y = y s . 



III. LOW- AND HIGH-ENERGY COMPONENTS 

Our constraints on the chiral pion-pion parameters involve integrals over absorptive 
parts. We have good experimental information in the range 600 MeV < y/s < 2 GeV but 
large uncertainties below this interval. As we cannot rely on experimental data in this low 
energy region we need some theoretical input. One way is to exploit the fact that the chiral 
amplitudes are meant to describe low energy scattering correctly. Therefore the low energy 
chiral absorptive parts can be used as approximations of the full absorptive parts. This idea 
was proposed ini and we apply it here. A more refined procedure is applied in! where the 
absorptive parts are computed from unitarized chiral S- and P-wave amplitudes. 
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In a first step the dispersion integral in (2.5) is split into a low energy integral Lj, from 
a = 4 to a = A 2 and a high energy integral, a > A 2 , H^. The representation (2.5) becomes 

Gi(x, y) = (l- 5 i2 )Gi(x , y ) + L t (x, y) + Hi(x, y). (3.1) 

This decomposition induces a crossing-symmetric decomposition of the amplitudes T 1 into 
low- and high-energy components. 

We assume now that A 2 can be chosen in such a way that the 2n-th order chiral absorptive 
parts A^(<7, r) approximate the full absorptive parts for 4 < a < A 2 up to 0(p 2( - n+1 ^) 
corrections. Furthermore we have already assumed in Section II that Gi is approximated 
by Gf in the neighbourhood of the symmetry point. Thus, if (x,y) is close to (x s ,y s ), our 
assumptions allow us to rewrite eq. (3.1) as follows: 

Gf(x, y) = (1 - 6 i2 )Gi(x , y ) + Lf(x, y) + H t (x, y) + 0(p 2 ^). (3.2) 

Hi is given by the integral in (2.5) extended over a > A 2 ; Lf is given by the same integral 
restricted to 4 < a < A 2 and with Bi replaced by Bf . 

A key observation is that Gf and Lf have the same discontinuities across their cuts for 
4 < a < A 2 . We show in the Appendix that, up to 6th order, this implies that the difference 
(Gf — L/f) has the following form: 

Gf(x, y) - Lf(x, y) = P^x, y) + H?(x, y). (3.3) 

Pi(x,y) is a polynomial of first degree and Hf is a high-energy component given by a 
dispersion integral starting at a = A 2 . Equation (3.2) becomes 

Pi(x, y) + H?(x, y) = (l- S^G^Xo, y ) + H t (x, y) + 0(p 2 ^). (3.4) 

The left-hand side is entirely determined by the chiral amplitudes whereas the right-hand 
side is determined, apart from Gi(xo,yo), by absorptive parts above A 2 . 

Those equations (2.3) and (2.4) which do not involve Gi(xo,yo) are now replaced by 
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(3.5) 



(P a + H*) {4) (x„y a ) = H 2 (x s ,y s ) + 0(p 6 ), 
— (P + tf *) (4) (x„ y s ) = ^(x s , y s ) + 0(p 6 ). 

(P 2 + tf 2 *) (6) (z„ = # 2 (x s , y s ) + 0(p 8 ), 
^(P i + H?)®(x s ,y s ) = -^H i (x s ,y s ) + 0(p 8 ), i = 0,l, (3.6) 
|- (P + tf*) (6) (x.,y.) = ^-H (x s ,y s ) + 0(p 8 ). 

The two equations (3.5) constrain the four 4th-order parameters and the four equa- 
tions (3.6) constrain the six 6th-order parameters. These equations clearly reduce to eqs (2.3) 
and (2.4) if A 2 = 4. Both sides of the equations depend on A 2 . However, the restrictions on 
the chiral parameters they imply should not depend on A 2 , up to higher-order corrections, 
as long as A 2 is in an interval where our assumption on the absorptive parts is valid. Thus 
we see that eqs (3.5) and (3.6) provide at the same time a tool for the determination of the 
chiral parameters and a verification of the validity of the chiral expansion. We elaborate on 
these two aspects in the following Sections. 



IV. INPUTS FOR THE PARAMETER CONSTRAINTS 

After these lengthy preliminaries we are ready for a feasibility test of our programme. 
To this end we need an Ansatz for the pion-pion absorptive parts and explicit forms for 
the chiral pion-pion amplitudes. We start with a definition of an Ansatz modelling the 
absorptive parts down to the elastic threshold. 

For 4 < s < 51 (280 MeV < \/s < 1 GeV) we retain only S- and P-wave contributions. 
A parametrisation proposed by Schenki is used for the corresponding phase shifts: 

1 



tan*/ (s ) = - 4))' (i-i) 3 i-ii [af - i( S - 4)11 

where I = for / = 0, 2 and I — 1 for I — 1, aj is a scattering length and 



(4.1) 



% = <t-^tl + («0 3 ' ( 4 - 2 ) 
10 



where bj is an effective range. In (4.1), si = M 2 = 30.4 and we adopt one of Schenk's choices 
for the remaining poles: s = 38 = (865 MeV) 2 , s 2 = -43 = -(920 MeV) 2 . Scattering 
lengths and effective ranges are chosen in accordance with the results of 4th-order standard 
chiral perturbation theoryi: 



a° = 0.20, a 2 = -0.042, a\ = 0.037, 



b° = 0.24, b 2 = -0.075, b\ = 0.005. 



(4.3) 



For s > 51 we adopt an Ansatz which has been used recently inB. In the interval 
51 < s < 110 (1 GeV < yfs < 1.5 GeV) we take only the /o contribution 



A%{s,t)=b*T h M h , 



■^M^M^Ml^)^-^ <4 ' 4) 
with T /o = 0.896 (= 125 MeV), M fo = 9.092 (= 1269 MeV). 

The high-energy region s > 110 is described in terms of Pomeron exchange and a degen- 
erate (p + /o) Regge 

trajectory^. If A^(s,t) desi gnates the isospin / t-channel absorptive 
part, the Pomeron exchange is specified by 



32ttM 2 Vx , 

with x = 10 GeV 2 , b = 10 GeV~ 2 , a' P = 0.4 GeV" 2 , A^ om = A^ m = 0§. The (p + f ) 
exchange is represented by the following absorptive parts: 

^l^t) = 4io( s ^ = 0.6sin(7ra(t)) T(l - a(t)) J (4.6) 

where a(t) = 1/2 + t/(2M 2 ) and Af} = 0. 

The question we are now asking is whether chiral perturbation theory can produce am- 
plitudes which are consistent with the absorptive parts defined in (4.1), (4.4), (4.5) and (4.6) 
in the sense that they satisfy the conditions (3.5) and (3.6) in a suitable range of values of 
A 2 . 

According to eqs ( |1.4| - |1.6D and (2.2) the totally symmetric one- and two-loop chiral 
amplitudes entering into our constraints are of the form 
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Gi(s,t,u) 



Qo(s,t,u) + Pa(s)f a (s) + U a (t)f a (t) + U a (u)f a (u)} 



a=l 
5 



Qx(s,t,u) + {3 [T a {s)f a (s) + T a (t)f a (t) + T a (u)f a \ 



u 



a=l 



+ 



ls-t 

5 



a=l ^ 



(W a (s)f a (s) - W a (t)f a (t)) + permutations 
1 



t — u 



(W a (t)f a (t)-W a (u)f a (u)) 



u — s 



(W a (u)f a (u)-W a (s)f a (s)) 



+ permutations | . 



(4.7) 



In these expressions Qo and Q\ are polynomials, Q2 is a constant, U a , T a and W a are 
polynomials obtained from those in (|1.6|) whereas the f a are the analytic functions appearing 
there. In standard chiral perturbation theory the Qi, U a , T a and W a depend linearly on the 
parameters b\, . . . , 6g. 



V. EVALUATING 4TH- AND 6TH-ORDER PARAMETERS 

The ingredients collected in the last Section allow us to compute the left- and right-hand 
sides of eqs (3.5) and (3.6). The 4th-order equations (3.5) involve only two parameters, 63 
and 64: 

3 



32tt 



A 



1 (41 

-63 + 36 4 + Cx(A 2 ) = H 2 (A 2 ) + 0(p 6 ), 



1 ' r 63 + 36 4 + C 2 (A 2 )l (4) =^(A 2 )+0( P 6 

J ox 



(5.1) 



-A 4 

7T 



The A 2 dependencies are indicated explicitly and every reference to the symmetry point 
is dropped. For instance H 2 (A 2 ) is the A 2 -dependent value of H2(x s ,y s ) and C*i(A 2 ) is the 
term of (P2 + H2)^(x 8 ,y 8 ) which is does not depend on the b{. 

For the time being it is convenient to work with the combinations (—63+364) and (63+364) 
rather than 63 and 6 4 . Figure 1 displays the values obtained form eqs (5.1) ignoring the 0(p 6 ) 
corrections. The variations of H 2 and d x H as functions of A 2 are compensated to a large 
extent by the A 2 -dependence of C\ and C 2 , but a sizeable A 2 -dependence remains. One may 
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ask if it could be compensated by A 2 -dependent 0(p 6 ) corrections. It is delicate to guess the 
order of magnitude of such corrections but we do not have to speculate at this point: as we 
shall see, the 0(p 6 ) corrections we are looking for are obtained from the 6th-order equations 
(3.6). 

Two of the equations (3.6) are corrected versions of the 4th-order eqs (5.1). They have 
the form: 



A 4 (-6 3 + 36 4 + Ci(A 2 )) (6) + A 6 [-4(6 5 - b 6 ) + 9l (b u b 2 , b 3 , h, A 2 ) (4) ] = ^H 2 (A 2 ) + 0(p 8 ), 



32tt 



(5.2) 



A 4 (63 + 36 4 + C 2 (A 2 )) (6) + A 6 [66 5 - 3b G + g 2 (b u b 2 , b 3 , 64, A 2 ) (4) ] = tt^(A 2 ) + 0(p 8 ) 



8Hn 



dx 



In these equations 

bf ) = bf ] + \ 2 5bi, i = 3, 4, cf ] = Cf + \ 2 5C k , k = 1,2, 

where 5bi and 5Ck are 6th-order contributions and 5Ck is known. 
In addition to (5.3) we have two true 6th-order equations: 



(5.3) 



A 4 Cf } (A 2 ) + A 6 [65 - 36 6 + g 3 (b u 62, 63, 64, A 2 ) (4) 
A 4 Ci 6) (A 2 ) + A 6 [3(6 5 + be) + <7 4 (&i, b 2 , 63, 64, A 2 )^ 



71 8Hq 2 



6 9j/ 



(A 2 ) + 0(p s ), 



2^(^+0^). 



(5.4) 



The fc in (5.2) and (5.4) are known linear and homogeneous combinations of 61, b 2 , b 3 and 
64 with A 2 -dependent coefficients. Fourth-order values have to be inserted. 

The constants C3 and C4 are of the form indicated in (5.3): they are sums of 4th-order 
terms and 6th-order corrections. We meet here a peculiarity of our programme. Although 
equations (5.4) do not appear at 4th-order level because they come from 0(p 6 ) terms in 
the expansions of the at the symmetry point, genuine 4th-order terms enter into these 
equations. This happens because each coefficient in the Taylor series of Hf has an expansion 
starting with a 4th-order term, in contrast with dPo/dy and dPi/dx which contain no such 
terms. In practice, the presence of 4th-order contributions in eqs (5.4) is harmless because 
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these contributions are small, especially for large enough A 2 , so that they have effectively 
the same order of magnitude as the true 6th-order terms. 

Eliminating the two 6th-order parameters 65 and b$ from eqs (5.2) and (5.4) gives two 
equations relating b\, b 2 , b 3 and 64: 



A 4 (-63 + 36 4 + C 5 (A 2 )) (6) + A 6 ^, b 2 , 63, h, A 2 ) 



(4) 



7T 

3 



dH, 



32if 2 (A 2 ) + ^(A 2 ) + 12^(A 2 ) 



dy 



dx 



+ 0(p 8 



(5.5) 



A 4 (63 + 36 4 + a(A 2 )) (6) + \ 6 g e (h, b 2 , b 3 , 64, A 



2\(4) 



71 

3 



-^o (a2) _^o (a2) _ 12 ^i (a2 , 
' ox 



+ 0(p 8 ). 



dx dy 

As the 4th-order constraints (5.1) do not provide well defined values for b^ and b^ we 



insert bf^ and b^' into #5 and g 6 . This amounts to a redefinition of the 0(p s ) corrections. 
We retain b^f 1 and b^ and drop the indices labelling orders in the chiral expansion from now 
on. To illustrate the quantitative content of eqs (5.5), we write them explicitly for A 2 = 10, 
without mention of the 0(p 8 ) corrections: 

- (1 - [0.14])6 3 + (3 + [0.24])& 4 + 4.17 • 10" 4 

-[6.12 • l(T 3 ]6i - [4.26 • 10~ 3 ]& 2 + [4.78 • 1(T 4 ] = 1.81 • 1(T 2 

(5.6) 

(1 - [0.12])& 3 + (3 + [0.10])6 4 - 3.18 • 10" 3 

+ [1.39 • 10" 3 ]&i - [1.29 • 1(T> 2 - [4.56 • 1(T 4 ] = 7.28 • 1(T 3 . 

Both sides of eqs (5.5) have been divided by A 4 . The square brackets indicate two-loop 
terms: numbers without square brackets are 4th-order terms. 

For A 2 = 10, the 4th-order version of eqs (5.6), i.e. eqs (5.1), is 



;( 6 ) 



- h + 36 4 - 1-52 • 10~ 3 = 1.54 • 10" 



(5.7) 



63 + 36 4 - 6.50 • 10" 4 = 1.15 • 10 



1-2 
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These equations are not identical to eqs (5.6) with the square brackets dropped. This is due 
to the 4th-order contributions of dHft jdy and dH*/dx, mentioned previously, in eqs (5.6). 

Equations (5.5) are treated as determining 63 and 64 in terms of 61 and 62, using the 
values of these parameters obtained from 1-loop standard perturbation theory^ 2 1 ! - 



6, 1 



96vr 2 
1 



13 

Mi " 3l 3 " I2Z4 + y 



-2/i + 3l 4 - ^ 



(5.8) 



24tt 2 

With the central valuesS'0 h = -1.7, T 3 = 2.9, f 4 = 4.6 we get 

6j = -7.7- 10" 2 , 6 2 = 7.1-10~ 2 . (5.9) 

Inserting these values into the A 2 = 10 equations (5.6) gives 

(-63 + 3& 4 ) (6) = 1-63 ■ 10~ 2 , (63 + 36 4 ) (6) = 1-12 - 10" 2 . (5.10) 

The corresponding 4th-order values are 1.70 ■ 10~ 2 and 1.21 • 10~ 2 , quite close to the 
6th-order values (5.10). The variation of (+63 + 364) at 4th and 6th order as functions of A 2 
obtained from eqs (5.1) and (5.5) is shown in Fig. 1. The striking feature is that (^63 + 364)^ 
is independent of A 2 within 4% for 4 < A 2 < 16. On the other hand, the difference between 
the 6th- and 4th-order values can be substantial, especially for (63 + 364) and A 2 between 4 
and 10. 

The values of 63 and 64 as obtained from (— 63 + 364) and (63 + 364) are displayed in Fig. 2: 
6 4 4 ^ and 6 4 6 ^ are practically constant, their difference due to a 6th-order correction is of the 
order of 4%. The parameter 63 is not as robust as 64. The way it is obtained amplifies the 
uncertainties. At 4th order, they must have the same magnitude as the computed values. 
The situation improves at 6th order and 63 becomes independent of A 2 between 4 and 16 
within 6%. 

Ignoring all uncertainties in the input quantities Hi and taking the variations with A 2 
for 4 < A 2 < 16 into account leads to the following 6th-order values of 63 and 64 : 

15 



b 3 = (-2.55 ±0.20) ■ 10 



-3 



6 4 = (4.55 ±0.15) ■ 10 



,-3 



(5.11) 



The only clear 4th-order result is 64 = (4.85 ± 0.15) • 10~ 3 . 

Once 63 and 64 are known, eqs (5.4) determine 65 and b$. With (5.9) and the central 
values in (5.11), this leads to depending weakly on A 2 , as shown in Fig. 3. In contrast, 
65 exhibits a strong A 2 -dependence: it decreases from 1.4 • 10~ 4 to 0.66 • 10~ 4 as A 2 varies 
between 4 and 16. Comparing Fig. 2 and Fig. 3, we observe that bf^ and b^ behave similarly, 
and so do b^ and 65 . This is possibly due to the similar roles which the pairs (63, 64) and 
(65,65) play m A pol (s;t,u) (eq. (|1.5|) ). Furthermore, the 4th-order absorptive parts do not 
depend on (63,64) and the 6th-order absorptive parts do not depend on (65, 6 6 ). One may 
expect that 65 would become stable at 8th order in the same way as 6 3 is stabilized at 6th 
order. 

Ignoring again the uncertainties in the Hi, Fig. 3 gives 



The values of 63, 64 and 6 6 indicated in (5.11) and (5.12) are compatible with the results 
obtained mi. Concerning 6 5 , no large instability shows up in@. This is probably related to 
the use of unitarized absorptive parts. In any case the order of magnitude we obtain is the 
same as in@. 



The last Section provides evidence that our programme is working. It produces three 
6th-order parameters which are stable in the sense that they depend weakly on A, the energy 
separating low and high energies as long as 2 < A < 4. The following comments summarise 
the salient features of our findings. 

1. The fact that 63, 6 4 and b§ are nearly independent of A 2 at 6th order is not the result 
of a weak A 2 -dependence of the individual terms in eqs. (5.2) and (5.4). In fact they vary 
strongly with A 2 : for instance H 2 (A 2 ) in eq. (5.2) increases from 1.01 ■ 10~ 2 to 1.69 ■ 10 -2 



6 6 = (9.2 ±0.3) • 10~ 5 . 



(5.12) 



VI. DISCUSSION 
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when A 2 goes from 4 to 16. Very delicate compensations are at work and are responsible for 
practically constant 6th-order parameters. 

2. The choice of the values (5.8) for b\ and 6 2 plays a role in the constancy of bf\ bf^ 
and bf\ Substantially different values would destroy the compensations mentioned under 
1. In fact, the constancy of bf^ and bf^ at the level displayed in Fig. 2 is certainly lost if 
(6i,6 2 ) is outside the square (-2 • 1CT 1 , 0) x (0, 2 • 1CT 1 ). 

3. One of our basic assumptions is that low-energy absorptive parts can be approximated 
by chiral absorptive parts. To check this we have computed the 4th- and 6th-order S- and P- 
wave absorptive parts. The 4th-order absorptive parts do not depend on the bf the central 
values (5.9) and (5.11) have been used in the 6th-order calculations. Figure 4 shows that in 
the case of Im /q(s) there is good agreement between the 6th-order chiral absorptive part 
and our Ansatz (4.1) if 4 < s < 18. There is also a spectacular improvement with respect 
to the 4th-order result. The agreement is not as good for Im fl(s) and Im /q(s). This 
does not affect our results because these absorptive parts are small at low energies and the 
contributions of Im /q (s) dominate. 

4. According to eqs (5.5), b 3 and 6 4 depend additively on Ci, C 2 , b 2 , H 2 , dH /dx : 
dHi/dx and dH /dy. This means that their values are dictated not only by the inputs 
H2, dH Q /dx, . . . but also by the structure of the chiral amplitudes. The inputs produce 
contributions given in Table 1 for four values of A 2 . This table also displays the relative sizes 
of the contributions coming from the three energy intervals (A, 1 GeV), (1 GeV, 1.5 GeV) 
and (1.5 GeV, 00) appearing in the input defined in Section IV. In this table, 63 [si, s 2 ] and 
M 5 !' s 2] are the contributions of the Hi to bf^ and bf^ coming from integrals restricted to 
Si < a < s 2 . 

We are working here with a crude Ansatz for the absorptive parts and we are primarily 
interested in estimates of the bi rather than precise determinations including error evalua- 
tions. Table 1 indicates how our input should be improved to favour a precise determination 
of 63 and 64. Although the major parts of 63 [A 2 , 00] and 64 [A 2 , 00] come from energies be- 
low 1 GeV, the energies above 1 GeV contribute to 64 [A 2 , 00] at the 10% level and there 
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is a strong cancellation at work for &3[A 2 ,oo]. We use a Regge Ansatz on (1.5 GeV, oo) 
and it would be difficult to estimate the errors associated with this. The cut at s — 110 
(= (1.5 GeV) 2 ) should be replaced by a cut at s — 200 (= (2 GeV) 2 ), experimental data 
used on (51, 200), the Regge Ansatz being restricted to (200, oo). This would approximately 
halve the Regge contribution. 

5. The implications of our findings for the 6th-order coupling constants are beyond the 
scope of the present work. We observe only that the 4th-order relation 

l 2 = ^.^.bf + l (6.1) 
6 

allow a determination of I2 (notice that I2 does not enter into our evaluation of b\ and 62, 
eqs (5.8)). This gives I2 = 5.4 which is below the central value 1% = 6.1 obtained from the 
analysis of A^ 4 decaysB A low value such as this has also been extracted from pion-pion 
scattering ini andEI. Still at 4th order we have 



b 



(4) = _J_ 
' 967T 2 



2/1 + h - ~ 



(6.2) 



With the value T\ = —1.7 used in eqs (5.8), and l 2 = 5.4, this gives 63 = — 1.6 • 10~ 3 . This is 
compatible with the range of variation of b^ displayed in Fig. 2 and confirms the credibility 
of our approach. 
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DETERMINATION OF LOW-ENERGY POLYNOMIALS 

We have to show that the terms Pj in eq. (3.3) are polynomials. We do this in detail for 



18 



The functions f a vanish at the origin and can be obtained from a once-subtracted dis- 
persion relation: 

s r°° dcr 1 



/«(*) 



-Im f a (a). 



(1) 



7r J4 a a — s 

By writing the right-hand side as the sum of a low-energy and a high-energy integral, f a is 
decomposed into a low- and a high-energy component. Introducing this decomposition into 
the expression (4.7) of G gives 

rA 2 da 



1 r l 

G%(s,t,u) = Q (s,t,u) + / 



a 



*t/ a (*) ull a (u) 
a — s a — t a — u 



Im f a (a)+H$(s,t,u). 

(2) 



The high-energy component ifg' is the sum of integrals in (0) extended to [A 2 , oo). The 
integrals in (Q) can be rewritten as 
1 r A2 da 



1 f da 

— / — [U a (s, a) + U a (t, a) + C/ a (u, a)] Im / a (a) 

71" J 4 a 

U a (a)lm f a (a) 



1 /-a 2 da 
+ - / — 

it J a a 



1 1 

+ 



a — s a — t a — it 



(3) 



where 



t/ a (s,a) 



sf/„(V) - aUJa) 



s — a 



(4) 



is a polynomial. On the other hand, (4.7) and (2.7) tell us that Bq = ^£?o, a with 

a 

BUs, t) = (s- t)(2s - 4 + t)U a (s)lm f a (s). (5) 



According to (2.5) and (2.6), we have Lq = J2 a Lo a , where 



£o,a(s>*,M) 



l r A2 



7T J4 



da 



So X ,a(^r). (6) 



(a - s)(a - t)(a - w) (a - s )(a - £ )(c - «o)_ 
It turns out that the difference between the second integral in (|3|) and Lq a is equal to 
this integral with (s,t,u) replaced by (so,t ,uo). Consequently (3.3) is true for i = with 

1 r A2 da 

P (s, t, u) = Q (s, t,u) V / — [U a (s, a) + U a (t, a) + U a (u, a)] Im f a (a) 

U a (a)lm f a (a). (7) 



1 f A \ 
+ - / da 

7T J 4 



1 1 

+ 



a — s a — t a — it 



19 



In terms of x and y, and at fixed (xo,yo), Po(x,y) is a polynomial of first degree. 
Similarly, 

1 f A2 da 

Pi(s,t,u) = Qx{s,t,u) V / — {3 [T a (s,a) +T a (t,a) +T a (u, a)] 

IT ~ J 4 O 



+ 



s-t 



(W a (s, a) — W a (t, a)) + permutations 



+ const. 



p 2 = Q2 _IW^< 1 



7T „ J a a Is — t 

1 



1 



u — s 



t — u 

(W a (u,a)-W a (s,a)) 



(W a (t,a)-W a (u,a)) 



permutations | 



The polynomials W a (s, a) and T a (s, a) are obtained in the same way as U a (s, a) 
It turns out that Pi(x,y) is linear in x, independent of y and P 2 is a constant. 
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FIGURES 

FIG. 1. The variations of the combinations (±63 + 364) as functions of A 2 . The 4th-order values 
(dashed lines) are obtained from eqs (5.1) and the 6th-order values (full lines) come from eqs (5.5) 
with 61 = -7.7 • 1(T 2 , 6 2 = 7.1 • 10- 2 . 

FIG. 2. The values of 63 and 64 at 4th order (dashed lines) and 6th order (full lines) resulting 
from Fig. 1. 

FIG. 3. The A 2 -dependent values of 65 and 6 6 obtained from eqs (5.4) with 61 = -7.7 • 1CT 2 , 
6 2 = 7.1 • 1(T 2 , 63 = -2.55 • 1(T 3 and 64 = 4.55 • 1(T 3 . 

FIG. 4. Comparison of the S- and P-wave absorptive parts defined by the Ansatz (4.1) (full 
lines) with the chiral 4th-order (dashed lines) and 6th-order (dashed-dotted lines) absorptive parts. 
The latter are obtained from eq. (1.6) with the same values of 61, 62, 63 and 64 as in Fig. 3. 
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TABLES 

TABLE I. The A 2 -dependent contributions to 63 and 64 at 6th order coming from the input 
quantities H2, dH^/dx, dH±/dx and dHo/dy according to eqs (5.5). The full contributions are 
given in the second column. The relative sizes of the contributions of the three energy intervals 
A 2 < s < 51, 51 < s < 110 and s > 110 appear in columns 3-5. 



A 2 



63 [A 2 , 



00 



h [A 2 , 51] 
b 3 [A 2 , 00] 



63 [51, 110] 
b 3 [A 2 , 00] 



63 [HO, 00] 
63 [A 2 , 00] 



12 
16 
A 2 



-5.94 • 10~ 3 
-5.78 • 10~ 3 
-7.60 • IO- 3 
-8.74 • 10~ 3 



64 [A 2 , 00] 



1.011 
1.004 
1.004 
1.007 



6 4 [A 2 , 51] 



64 [A 2 , 



00 



-0.060 
-0.060 
-0.049 
-0.046 



64 [51, 110] 
64 [A 2 , 00] 



0.049 
0.055 
0.045 
0.038 



64 [no, 



DC 



& 4 [A 2 



00 



12 
16 



4.94 • 10~ 3 
4.49 • 10~ 3 
4.02 • 10~ 3 
3.37- 10~ 3 



0.905 
0.900 
0.894 
0.888 



0.017 
0.018 
0.019 
0.019 



0.078 
0.082 
0.088 
0.093 
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